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Abstract. We consider generation and evolution of small-scale magnetic fields in neutron stars. These fields can be generated 
by small-scale turbulent dynamo action soon after the collapse when the proto-neutron star is subject to convective and neutron 
finger instabilities. After instabilities stop, small-scale fields should be frozen into the crust that forms initially at high density 
~ 10^'' g/cm^ and then spreads to the surface. Because of high crustal conductivity, magnetic fields with the lengthscale ~ 1 — 3 
km can survive in the crust as long as 10-100 Myr and form a sunspot-like structure at the surface of radiopulsars. 
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1. Introduction 

Our knowledge of the magnetic field strength in neutron stars 
comes mainly from radiopulsars with measured spin-down 
rates. With the assumption that the spin-down torque on the 
pulsar is determined by its magnetodipole radiation (Ostriker 
& Gunn 1969), the spin period, P, and spin-down rate, P, are 
related to the field strength at the magnetic pole, Bd, by 



PP 



(1) 



where / is the moment of inertia and R is the radius (we as- 
sume that the magnetic and rotation axes are perpendicular). 
The magnetic fields inferred from the spin-down data range 
from ~ 5 X lO^'^ to ^ lO^G but, most likely, these fields charac- 
terize the global magnetic configuration of neutron stars rather 
than a fine magnetic structure near the stellar surface. 

Measurements of the spin-down rate are not, however, the 
only way to obtain information about the neutron star mag- 
netic fields. Recent observations of the X-ray spectra features 
of some pulsars provide one more opportunity to look into the 
magnetic field near the surface. Absorbtion features in the spec- 
trum of the isolated pulsar IE 1207.4-5209 were associated by 
Sanwal et al. (2002) with atomic transitions of once-ionized 
helium in the neutron star atmosphere with a strong "surface" 
magnetic field, Bg ^ 1.5 x lO^^G. An estimate of the spin- 
down rate of this pulsar provides the strength of the "dipole" 
magnetic field, B^ (2-4) x IQi^G (Pavlov et al. 2002) that is 
typical for radio pulsars of that age 0.2—1.6 Myr). Becker et 
al. (2002) found an emission line in the X-ray spectrum of PSR 
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B 182 1-24 that could be interpreted as cyclotron emission from 
the corona above the pulsar's polar cap. This emission Une is 
likely formed in a magnetic field Bs ^ 3 x lO^^G which is ap- 
proximately two orders of magnitude stronger than the "dipole" 
field inferred from P and P. Haberl et al. (2003) interpreted a 
broad absorption feature in the spectrum of the isolated neu- 
tron star RBS1223 as a cyclotron line produced by protons in 
the magnetic field ~ (2 — 6) x 10^"^ G. These measurements 
provide strong evidence that the local magnetic fields on the 
neutron star surface can exceed the conventional "dipole" field. 

Another piece of evidence on the distinction between the 
"dipole" and "surface" magnetic fields comes from the data on 
emission of radio pulsars. Recently, Gil & Mitra (2001) and 
Gil & Melikidze (2002) argued that the formation of a vacuum 
gap in radio pulsars is possible if the actual surface magnetic 
field near the polar cap is very strong, Bg ^ lO^'^G, irrespec- 
tive of the field measured from the P-P data. Also, radio emis- 
sion from the recently discovered pulsar PSR J2 144-3933 with 
the longest period 8.5 s, which lies extremely far beyond the 
conventional death line, can be understood if this pulsar has 
a strong surface magnetic field of complex geometry (Gil & 
Mitra 2001). Radio emission of many other pulsars which lie 
in the pulsar graveyard and should be radio silent can be ex- 
plained if one adopts the model with a strong and complex 
surface field with a small curvature of the field lines (< 10^ 
cm). This model is consistent with the conclusion of Arons & 
Scharlemann (1979) and Arons (1993) that pulsars with very 
long periods (> 5 s) require a more complex field configura- 
tion than a dipole if pair creation is essential for the mecha- 
nism of radio emission. Analysing the phenomenon of drift- 
ing subpulses observed in many pulsars, Gil & Sendyk (2000) 
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found that their behaviour is consistent with the vacuum gap 
maintained by a strong sunspot-like magnetic field. Following 
this idea Gil, Melikidze & Mitra (2002 b) found that in the fa- 
mous case of PSR 0943+10 (Deshpande & Rankin 1999, 2001) 
the "surface" field should be approximately 7 times stronger 
than the "dipole" component inferred from the spin-down rate. 
A sunspot-like configuration of the surface magnetic field is 
also suggested by the spin-down index in some pulsars (Cheng 
& Ruderman 1993, Ruderman, Zhu & Cheng 1998). Also, 
Cheng & Zhang (1999), analysing the X-ray emission from 
the polar regions of the rotation-powered pulsars, argued that 
Bg ^ lO'^'^G and the characteristic curvature of the field fines 
is around lO^cm resembling the sunspot-like structure. 

Rapidly growing amount of evidence on the distinction be- 
tween the local field strength at the stellar surface and the 
global "dipole" field suggests that this can be a general phe- 
nomenon in neutron stars. In this paper, we propose a scenario 
of the formation of such complex magnetic configurations with 
the strength of a small scale "surface" field in excess of the 
"dipole" component. It is generally accepted that neutron stars 
are subject to hydrodynamic instabilities soon after their birth 
in the core collapse (see, e.g., Epstein 1979, Livio, Buchler 
& Colgate 1980, Burrows & Lattimer 1986). The convective 
stage lasts about 30-40 s (Miralles, Pons & Urpin 2000, 2002) 
and, under certain conditions, turbulent motions can amplify 
the magnetic field via dynamo action (Thompson & Duncan 
1993, Xu & Busse 2001). The generated magnetic field must be 
frozen into the crust that is formed in the course of neutron star 
cooling. Because of high crustal conductivity, magnetic fields 
with a relatively small lengthscale, ^ 10^ cm, can survive dur- 
ing the active lifetime of radio pulsars. 

The paper is organized as follows. In Section 2, we consider 
convection and associated dynamo action in protoneutron stars, 
and estimate the field that can be generated during this stage. 
We discuss the crust formation and small scale crustal magnetic 
structures in Section 3. Our results are briefly summarized in 
Section 4. 

2. Convection in protoneutron stars 

Neutron stars are formed in a supernova explosion associated to 
the gravitational collapse of massive stars. The explosion lasts 
~ 1 s and is followed by the core bounce and generation of a 
strong shock which heats the protoneutron star to a very high 
temperature ^ 10^^ K (Burrows & Lattimer 1986, Burrows & 
Fryxell 1992, Rampp & Janka 2000). Hydrostatic equilibrium 
settles down on a very short timescale ^ 10""^ — 10~^ s but, 
even when the equilibrium is reached, the surface layers are 
very extended, thus the radius of a protostar is ~ 50 — 100 km 
instead of the canonical 10-15 km. However, compression of 
the exterior zones is quite fast and, after a few seconds, the star 
is compressed to the canonical radius. 

Hydrodynamic instabiUties in a newly born neutron star 
are driven by both the lepton gradient (Epstein 1979), and the 
development of negative entropy gradients, which are com- 
mon in many simulations of supernovae explosions (Bruenn 
& Mezzacappa 1994, 1995, Rampp & Janka 2000) and evolu- 
tionary models of protoneutron stars (KeU & Janka 1995, Keil, 



Janka & MuUer 1996, Pons et al. 1999). Likely, both, the con- 
vective and neutron finger instabiUties, can arise in protoneu- 
tron stars (Miralles, Pons & Urpin 2000) with the neutron finger 
unstable region typically surrounding the convective region. 
Initially, only the surface layer, containing around 0.5Mq, is 
unstable but the bottom of the unstable region spreads down to 
the center. Approximately in 10 s, the whole star is subject to 
instability with convection, operating in the central region of 
the enclosed mass ^ 0.5Mq, and the neutron finger instability 
dominating in the rest of the volume. After ~ 20 s, the unstable 
region begins to shrink to the center and, at t 40 s, both the 
temperature and lepton gradients become too smooth to main- 
tain instabilities. Both, the Rayleigh and Grashoff numbers are 
typically large in the unstable regions, and instabiUties Ukely 
do operate in a turbulent regime (Thompson & Duncan 1993). 

During the hydrodynamic unstable phase, a protoneutron 
star is opaque to neutrino, and the turbulent velocity can be 
estimated by the standard mixing-length approximation (see, 
e.g., Schwarzschild 1958). The largest unstable lengthscale is 
of the order of the pressure lengthscale, L, and the turbulent 
velocity in this scale, vl, can be estimated as 



L 

TL ' 



(2) 



where tl is the instabiUty growth time (which is of the order 
of the turnover time in a scale L). In the convectively unstable 
region, we have 



1 1 1 JAvri 
-2 ~ = i^gP^f^: 



(3) 



where Tc is the convection growth time, g is the gravity, AVT 
is the difference between the actual and adiabatic temperature 
gradients, and (3 is the thermal expansion coefficient. The tur- 
bulent velocity, vl, varies in time since AVT progressively 
reduces due to the neutron star cooling. Convection is a dy- 
namical instability and grows on a short timescale (^ 0.1 — 1 
ms when convection is most efficient). Assuming L ^ 1 — 3 km 
and using the calculations of Tc (Miralles, Pons & Urpin 2000), 
we can estimate vl ^ 10^ cm/s (except the late unstable phase 
when convection is almost exhausted). 

In the region of a "doubly diffusive" instability often re- 
ferred to as "neutron fingers", a destabilizing influence of the 
lepton number gradient usually dominates the effect of AVT. 
This instability is the astrophysical analog of salt fingers that 
exist in terrestrial oceans. Physically, a fluid element perturbed 
downward in the proto-neutron star can thermally equilibrate 
more rapidly with the background but find itself lepton-poorer 
and denser and, therefore, subject to a downward force that 
would amplify perturbations. This instability is typically more 
efficient in the region above the convective zone, involving a 
larger portion of the stellar material. With the same reasoning, 
we can estimate vlih this region substituting 



1 



i- = \9S\VY\ 



(4) 



into equation (2); r„/ is the neutron finger instabiUty growth 
time, S is the chemical expansion coefficient and Y = (ue + 
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nv)/n is the lepton fraction with Ug, n^, and n is the number 
density of electrons, neutrinos, and baryons, respectively. The 
neutron finger instability is typically slower than convection, 
T„/ ~ 30 — 100 ms (except the very early and very late phases). 
Being slower, this instability can nevertheless exist in a more 
extended region. Using the calculations of t„/ (Miralles, Pons 
& Urpin 2000), we estimate ~ (1 - 3) x 10'' cm/s in the 
neutron finger unstable region. 

Likely, protoneutron stars rotate relatively rapidly (Zwerger 
& Miiller 1997, Rampp, Miiller & Ruffert 1998) and, hence, 
can be subject to the turbulent dynamo action. The initial spins 
of pulsars are not well constrained by observations but, most 
hkely, they lie around ~ 100 ms (Narayan 1987). The influence 
of rotation on turbulence is characterized by the Rossby num- 
ber, Ro= P/tl where P is the spin period. Since there are two 
unstable regions inside the proto-neutron star with substantially 
different properties, the Rossby number can differ much in the 
convective and neutron finger unstable regions. In the convec- 
tive zone, we have Ro~ 100, and the influence of rotation 
on turbulence is probably negUgible. On the contrary, in the 
neutron fingers unstable region, Ro^ 1 and turbulence can be 
strongly modified by the Coriohs force. Therefore, the neutron 
finger unstable region seems to be better suited for the mean- 
field dynamo action (see Bonanno, Rezzolla & Urpin 2003). 
Note that this region is typically above the convectively unsta- 
ble region, therefore the mean-field dynamo operates mainly in 
the surface layers. 

In both unstable regions, however, turbulent motions can 
generate turbulent magnetic fields by small-scale dynamo ac- 
tion. The electrical conductivity of hot nuclear matter, a, is rel- 
atively high, « 1.5 X lO^^T^o^ s-i where Tio = T/IQ^^K 
(Baym, Pethick & Sutherland 1971), and therefore the char- 
acteristic timescale of ohmic dissipation is very long (^ 10^ 
yrs for the field with the lengthscale L ~ 1 — 3 km). The 
magnetic diffusivity, rj = (? jAnxa, is small compared to vis- 
cosity, V, therefore the dynamo operates in the regime of large 
magnetic Prandtl numbers, Pr=i^/?7. There are two factors that 
cause viscous stresses in protoneutron stars: neutrino transport 
and neutron scattering (Thompson & Duncan 1993). Neutrino- 
induced viscosity dominates on scales larger than the neutrino 
mean-free path and can be very efficient. However, even this 
large viscosity cannot prevent instability on scales compara- 
ble to the pressure lengthscale (Miralles, Pons & Urpin 2000). 
The viscosity caused by neutron scattering operates on scales 
shorter than the neutrino mean-free path and is much smaller, 
but the magnetic Prandtl number is large even in this case, 

Pr = ^ « 2.6 X lOVii^Tio^ » 1, (5) 

where p\i = p/lO^^g/cm^ and p is the density; we use the an- 
alytical fit for the coefficient of viscosity of hot nuclear matter 
obtained by Cutler, Lindblom & Splinter (1990). 

There are two qualitatively different phases to a small-scale 
dynamo action: the kinematic phase when the field does not 
provide a noticeable influence on turbulent motions, and the dy- 
namical phase when velocity is affected by the Lorentz force. 
In the kinematic regime, a weak seed magnetic field grows in 
strength exponentially in time, while the characteristic length- 
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scale of a field, Ib, decreases exponentially (Kazantsev 1968, 
Kraichnan 1976, Kulsrud & Anderson 1992). In protoneutron 
stars, instabilities generate primary turbulent motions on scales 
~ L, and these motions amplify turbulent magnetic fields of 
the same scale. However, the energy of such magnetic fluctua- 
tions is transfered to the small scales after a few eddy-turnover 
times. For Pr:^ 1, the magnetic energy is first transfered to 
scales shorter than the viscous dissipative scale, £ Re-3/4i 
(Re is the Reynolds number at the lengthscale L), but greater 
than the magnetic dissipative lengthscale, in ~ £i/Pr~^/^ (see, 
e.g., Schekochihin, Boldyrev & Kulsrud 2002). Bofli lengfli- 
scales, £i, and are very small in protoneutron stars, 

4 « 1.9 X W-%'/'lI/%-,'/'pT' cm, (6) 

i, « 1.2 X 10-%'^'lI/X^'pT cm, (7) 

where vs = Ui/lO^cm/s, and L5 = L/lO^cm. Turbulent mo- 
tions with the lengthscale shorter than £1, are suppressed by vis- 
cosity, but magnetic fluctuations can exist if their lengthscale is 
larger than Probably, turbulence is well developed in both 
unstable zones, and fluctuations spread through a wide range 
of scales (see, e.g., Thompson & Duncan 1993). The kinematic 
growth phase is terminated when the Lorentz force starts mod- 
ifying the convective motions. Since the characteristic growth 
time of small-scale fields, ^ L/vL,is much shorter than the du- 
ration of the unstable phase, the small-scale dynamo operates 
likely in a nonUnear regime in both unstable regions. 

In the non-linear regime, turbulence drives kinetic and mag- 
netic energy cascades which are quasi-steady in protoneutron 
stars. The nature of these cascades in MHD turbulence has 
been a matter of debate for many years. Iroshnikov (1963) and 
Kraichnan (1965) suggested that in the inertial range kinetic 
and magnetic power spectra are given by 

Ek{£)(xEm{£)(x£^/^ (8) 

where i is the lengthscale of fluctuations. On the other hand, so- 
lar wind data suggest a spectrum cx £^^^ as in the Kolmogorov 
theory. Numerical simulations (see, e.g., Kida, Yanase & 
Mizushima 1991, Haugen, Brandenburg & Dobler 2003) indi- 
cate that the magnetic energy spectrum probably does not show 
a power law behaviour. It seems that spectra can be even shal- 
lower than ^^/^ in some models. For any more or less plausi- 
ble spectra, however, the amplitude of magnetic fluctuations in- 
creases with i and reaches its maximum on scales comparable 
to the main lengthscale of turbulence, L. In a saturated phase, 
small-scale magnetic fields are approximately in equipartition 
with velocity fluctuations, 

Bi ~ B,q{e) = s/^pvi, (9) 

where and vn are the ampUtudes of magnetic and veloc- 
ity fluctuations with the lengthscale £. Equation (9) yields only 
the order of magnitude estimate, since the magnetic and ki- 
netic energy have generally different spectra (Kida, Yanasa & 
Mizushima 1991). 

Turbulence is non- stationary in both unstable zones of pro- 
toneutron stars. It arises very rapidly soon after the coUapse, 
reaches some quasi-steady regime, and then goes down when 
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the temperature and lepton gradients are smoothed (after ~ 
30 — 40 s). The timescale required for fluid to make one turn 
in a turbulent cell with the lengthscale I can be estimated as 
« T^ljvi. This timescale varies with time, but is typically 
much shorter than the characteristic cooling timescale, TcooU 
except the very late phase when gradients are smoothed and in- 
stabilities are less efficient. Therefore, turbulence can be treated 
in a quasi-steady approximation during the almost whole un- 
stable phase. Using equation (9), we can estimate the maxi- 
mum field generated in the lengthscale I ^ L during the quasi- 
steady regime as Bj^c ~ 10^^ G in the convective zone, and 
^Lnf ^ (1 — 3) X 10^^ G in the neutron finger unstable zone. 
However, the temperature and lepton number gradients are pro- 
gressively reduced as the protoneutron star cools down and, 
therefore, the turbulent velocity decreases as well. As the re- 
sult, the strength of small-scale magnetic fields generated by 
turbulence also decreases compared to the maximum value, 
but estimate (9) is still valid until the quasi-steady condition 
Tcooi ^ T^, is fulfilled. We assume that this condition breaks 
down at f = fa when Tiify becomes comparable to the cooling 
timescale: Ti{ta) ~ Tcooi{ta)- Then, the turbulent velocity at 
t = ta is given by 



ve{ta) 



'Tcool {taY 



(10) 



We assume that the final strength of the magnetic field gener- 
ated by the small-scale dynamo, B^, is of the order of B^q {£) at 
t = ta- Perhaps this estimate provides only the upper hmit on 
a field strength since turbulent motions can influence magnetic 
fluctuations even at i > ta- However, the decay of turbulence 
is rather fast, and the kinetic energy of turbulence rapidly be- 
comes smaUer than the magnetic energy. Therefore, we have 
for the final strength of the generated magnetic field 



Bir 

where g is a numerical factor, q 
(10) for V(,{ta), we obtain 



(11) 

1. Substituting expression 



(12) 



The final strength of the generated small-scale field turns out 
to be the same for both unstable zones and decreases with de- 
creasing lengthscale. For the largest turbulent scale, £ ^ L ^ 
1 — 3 km, estimate (12) yields Bl ~ 3 x IQ^^q G if TcooI is 
of the order of a few seconds. Note that our estimate is in con- 
trast to that of Thompson & Duncan (1993) who assumed that 
the strength of the generated field is approximately given by its 
maximum value. 

3. Crust formation and crustal small-scale 
magnetic fields 

It is difficult to predict final disposition of the magnetic field 
after the turbulent motions inside the protoneutron star stop. 
Most likely, there exists a wide spectrum of magnetic fluctu- 
ations with the lengthscale shorter than L and with the field 
strength given by equation (12). When convection is exhausted. 



these small-scale fluctuations evolve under the influence of 
the ohmic dissipation and buoyancy. After a magnetohydrody- 
namic quasi-equilibrium is established, one expects that mag- 
netic loops densely fill the volume and surface of the star, al- 
though the field strength may vary considerably, depending on 
the degree of intermittency of the generated small-scale field. 
In liquid nuclear matter, the ohmic decay timescale of fluctua- 
tions with the lengthscale i is 



ro(^)~^ = 2.8xloVo'^^ 



yrs, 



(13) 



and dissipation is important only for fluctuations with very 
small £. 

Likely, the crust formation provides the most important in- 
fluence on the evolution of small-scale magnetic fields at this 
stage. Approximately at the age ^ 40 — 50 s (soon after con- 
vection stops), the neutron star cools down to the internal tem- 
perature ~ (1 - 3) X 10i°K (Pons et al. 1999). At such tem- 
perature, neutrons and protons can form nuclei and clusters in 
the nuclear matter with the density ^ 10^'* g/cm^. Wlien the 
neutron star cools down to a lower temperature, nuclei can be 
formed at p < 10^** g/cm^ as well. Note that nuclear composi- 
tion at high density depends generally on the pre-history of the 
neutron star. For instance, the composition of "ground state" 
matter (Negele & Vautherin 1973) differs noticeably from that 
of "accreted" matter (Haensel & Zdunik 1990). The Coulomb 
interaction of nuclei leads to the crystal formation even at a 
relatively high temperature. Crystallization occurs when the 
ion coupling parameter F = Z'^e^ /{akBT) reaches the criti- 
cal value r = r„ « 170 (Slattery, Doolen & De Witt 1980); 
a = the mean inter-ion distance, and Z are 

the number density and charge number of ions, respectively; T 
is the temperature, and fcs is the Boltzmann constant. Then, the 
crystaUization temperature is 



akB^m 



1.3xlO^Z^/3^-V3pV3lI2 K 



(14) 



where /ig is the number of baryons per one electron, and 
pi2 = p/lQ-^^ g/cm^. For the density p ^ 10^^ - 10^** g/cm^, 
the crystallization temperature is of the order of 10^" K (Baiko 
& Yakovlev 1996). Therefore, the crust formation starts almost 
immediately after the end of the convective phase, and the mag- 
netic fields generated by convective motions should be frozen 
into the crust. Solidification proceeds rather rapidly, and the 
outer boundary of the crust reaches the density ^ 10^° glcw? 
after ~ 1 day in the case of standard cooling. 

In crystal layers, the magnetic field evolution is mainly de- 
termined by ohmic dissipation. The crust electric conductivity 
(Tc, can be expressed in terms of the electron relaxation time, 
Te (see, e.g., Baiko & Yakovlev 1996), 



1.5 X IQ^^x^ 



(15) 



where x = pp/meC, (3 = vp/c; pp and are the electron 
Fermi momentum and velocity, respectively. In a high density 
region, electrons are ultrarelativistic and /3 w 1. Using expres- 
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sion (15), we can estimate the ohinic decay timescale of mag- 
netic fluctuations with the lengthscale ^ in the crust, 

TocW = « 0.7 X 10^x^(3 (y^Ss^) yrs, (16) 

where £5 = ^/lO^cm. This timescale is short for short length- 
scale fluctuations, but roc can be much longer than the cooUng 
timescale for fluctuations with relatively large L For instance, 
the decay of fluctuations with £ ~ 10^ — 10^ cm proceeds 
certainly on a timescale longer than the time required for the 
formation of a well developed crust. Therefore, turbulent mag- 
netic fields with relatively large lengthscales are likely frozen 
into the crystaUized matter soon after convection stops. 

Further evolution of turbulent magnetic fields deposed in 
the neutron star crust depends on conductive properties of the 
crustal matter and the cooUng scenario. The behaviour of such 
fields is governed by the standard induction equation where 
only the ohmic dissipation is included and is qualitatively sim- 
ilar to the behaviour of a large scale crustal magnetic field con- 
sidered in detail by Urpin & Konenkov (1997). We refer the 
results of this paper in what follows. We consider the evolution 
of turbulent fields for the neutron star model with the standard 
cooUng since this cooUng scenario can better account for the 
available observational data. Note that models with accelerated 
cooUng always lead to a slower decay of the magnetic field 
(Urpin & Van Riper 1993) and, as a result the small-scale mag- 
netic structure can survive for a longer time in such models. 
The crustal conductivity is determined by scattering of elec- 
trons on phonons and impurities. Scattering on phonons dom- 
inates the conductivity while the neutron star is relatively hot 
whereas impurities give the main contribution at a low crustal 
temperature (see, e.g., Yakovlev & Urpin 1980). The effect of 
impurities on the conductivity is characterized by the impurity 
parameter Q, 

Q=-y^n'{Z-Z'f, (17) 

n' 

where n' is the number density of an interloperspecies of 
charge Z' , n and Z are the number density and charge of 
the dominant ion species; summation is over all species. Most 
likely, Q ranges from 0.001 to 0.1 within the crust. 

In Fig.l, we plot the time dependence of the surface 
strength of a small-scale magnetic field for a neutron star model 
with mass M = 1.4Af[0] and with the equation of state of 
Pandharipande and Smith (see, e.g., Pandharipande, Pines 8c 
Smith 1976). The radius of this model is i? = l-'j.QS km. For 
the sake of simpUcity, we assume in calculations that the de- 
pendence of a field on polar and azimuthal coordinates is si- 
nusoidal with the wavelength equal to the main scale of turbu- 
lence, L, and use the so called local approximation in these two 
directions. This is a sufficiently good approximation since the 
radius of the star is substantially larger than L and the thick- 
ness of the crust. We also assume that the initial radial depth 
of a smaU-scale field is equal to the main scale of turbulence, 
L, and varies within the range 1-2 km. This range of depth 
corresponds to the density from w 5 x 10^^ g/cm^ to « 10^^ 
g/cm^ for the PS-model. The initial radial dependence of the 



T ' r 




Ig t, years 

Fig. 1. The evolution of the surface magnetic field strength B 
for different initial lengthscales L. The decay curves are shown 
for 1/ = 2 km (curve 1), 1.2 km (curves 2 and 4), and 1.0 km 
(curve 3). The initial field is 3 x lO^^ q ^j. 

curves 1, 2 and 

3, and 10^^ G for the curve 4. The crustal impurity parameter 
is Q = 0.01. 

field is chosen in accordance with the model proposed by Urpin 
& Konenkov (1997). Note that main conclusions of our paper 
are not sensitive to a particular choice of this dependence but 
are flexible to the value of L. 

The magnetic field evolution is shown for three initial 
depths of the field, L — 1 km (this depth corresponds to 
the density p « 5 x 10^^ g/cnr^ in the crust), L = 1.2 km 
ip PS 10" g/cm^), and L = 2km (p ss lO" km). For the 
purpose of illustration, the decay in the case L = 1.2 km is 
shown for two initial field strength, 3 x 10^^ G and 10^^ G 
(curves 2 and 4, respectively). During the initial stage (t < 10^ 
yrs), the crustal conductivity is determined by scattering of 
electrons on phonons and is relatively low for all calculated 
models. Therefore, the field decay can be essential during this 
stage, and the surface field strength can weaken by a factor 
« 2 — 7, depending on the lengthscale of the initial turbulent 
field. Obviously, a decrease during the initiale stage is smaller 
for the field with larger L. For instance, a smaU-scale field with 
L = 1 km is reduced by a factor rs 7 after 10^ yrs whereas the 
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field with L = 2 km decreases less than twice after the same 
time. Note that for L > 2 km the decrease during the initial 
stage is practically negligible, and the field at f ^ 10^ yrs is 
only a bit lower than the initial field. 

After ^ 10^ yrs, the dominant conductivity mechanism 
changes from electron-phonon to electron-impurity scattering. 
As a result, the conductivity increases and the rate of field 
decay slows down. We calculate the field decay during this 
stage using an intermediate value of the impurity parameter, 
Q = 0.01 that characterizes not a very poluted crust. For 
larger Q, the decay is more rapid (see, for comparison, Urpin & 
Konenkov 1997). The evolution of small-scale fields shows the 
presence of flat segments of the decay curves att~0.1 — 100 
Myr like those of a large-scale field. The length of a plateau de- 
pends on the impurity parameter, Q. The lower that Q gets, the 
longer the plateau on the corresponding decay curve is. During 
the impurity dominating stage, the decay turns out to be ex- 
tremely slow, and the characteristic decay time can be as long 
as 10- 100 Myr. 

These simple model calculations show very clearly that tur- 
bulent magnetic fields with the lengthscale of the order of the 
turbulence main lenghtscale in proto-neutron stars (~ 1 — 3 km) 
can survive in the crust during a very long time ^ 10 — 100 Myr 
that is generally comparable to the active lifetime of radiopul- 
sars. The field strength in such magnetic spots on the surface 
can reach 5 x 10^^ — 2 x 10^^ G depending on the radius of a 
spot and can be larger than (or comparable to) the strength of 
the dipole field. 

4. Conclusion 

We have considered the formation and evolution of small-scale 
magnetic structures in the neutron star surface layers. A wide 
spectrum of these structures can be generated by the small- 
scale turbulent dynamo action during the unstable phase that 
lasts ^ 30 — 40 s after the neutron star birth. There are two sub- 
stantially different unstable regions in the proto-neutron star, 
with the convective instability active in the inner region and 
the neutron-finger instability more efficient in the outer re- 
gion. Generally, the small-scale dynamo can generate small- 
scale magnetic structures in both unstable zones. Due to high 
conductivity of nuclear matter, the lengthscale of generated 
magnetic structures spreads from the main scale of turbulence, 
L ~ 1 — 3 km (comparable to the pressure scale height), to 
extremely short lengthscales determined by ohmic dissipation. 
After instabilities stop, structures with short lengthscales decay 
on a short timescale oc £^ due to finite electrical conductivity 
whereas fields with larger lengthscales can survive for a longer 
time. The crust formation that starts almost inmiediately after 
instabilities stop may have the decisive influence on the evolu- 
tion of such magnetic fields. These magnetic structures can be 
frozen into the crystallized matter and then evolve in the crustal 
layers. The neutron star cooling increases the conductivity of 
the crust and, as a result, small-scale magnetic structures decay 
extremely slowly. Our calculations show that structures with 
the lengthscale L ~ 1 — 3 km can survive as long as 10-100 
Myr that is basically comparable to the active life-time of ra- 
diopulsars. 



The strength of the magnetic field generated by the small- 
scale dynamo action is approximately determined by equipar- 
tition and decreases with the decreasing lengthscale. For mag- 
netic structures with L ~ 1 — 3 km, the magnetic field can be 
as strong as 5 x 10^^ — 2 x 10^^ G even for radiopulsars 
as old as ^ 10 — 100 Myr. Note that the decay of small-scale 
magnetic fields is qualitatively similar to that of the large-scale 
crustal field considered by Urpin & Konenkov (1997). The ori- 
gin of a large-scale field in neutron stars is still debatable, but 
it is possible that this field was generated by some mechanism 
in the layer with p < 2 x 10^^ g/cm^ that corresponds to the 
crust (see, e.g., Bonanno, Rezzolla & Urpin 2003). Then, the 
thickness of a layer occupied by the large-scale field can gen- 
erally differ from L. If the initial depth of a large-scale field is 
smaller than L, then the small-scale field decreases slower, and 
we observe a radiopulsar with magnetic "spots" where the field 
is stronger than the dipole field inferred from the spin-down 
rate. On the contrary, if the initial depth of a large-scale field 
is larger than L, then small-scale structures decay faster than 
the dipole field, and the resulting magnetic field becomes more 
regular with the age. 

Small-scale magnetic structures with L ~ 1 — 3 km and 
S ^ 5 X 10^^ - 2 X 10^^ G at the surface may have an impor- 
tant influence on many properties of radiopulsars. As already 
mentioned in the Introduction, this range of the surface mag- 
netic field is favourable for the inner vacuum gap formation 
in pulsars (Gil & Melikidze 2002). The sparking discharge of 
this gap produces filaments of electron-positron plasma, whose 
presence seems absolutely necessary for generation of coher- 
ent pulsar radio emission. Moreover, the E x B drift of spark 
plasma should be manifested as the observed subpulse drift, 
provided that the actual surface field has some degree of ax- 
ial symmetry, with a tendency to converge at the local pole 
(see Gil, Melikidze & Geppert, 2003 for review). These authors 
assumed that one small-scale structure (a "spot") coincides to 
some extent with the canonical (dipolar) polar cap in the sence 
that the magnetic field lines form a complex magnetic configu- 
ration near the pole but connect smoothly with a subset of open 
dipolar field lines at a larger radius (Gil, Melikidze & Mitra 
2002a). This means that the actual polar cap is defined by those 
non-dipolar field lines which penetrate the light cyUnder. 

A study of small-scale magnetic structures can also pro- 
vide information regarding the very early evolutionary stage 
of neutron stars since the observed fields were frozen into the 
crust soon after the neutron star birth. For example, the pres- 
ence of small-scale magnetic structures in radiopulsars can be a 
good evidence that proto-neutron stars pass the turbulent stage 
in their evolution. 
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